Geometric phase in the Kitaev honeycomb model and scahng behavior at critical 
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In this paper a geometric phase of the Kitaev honeycomb model is derived and proposed to 
characterize the topological quantum phase transition. The simultaneous rotation of two spins is 
crucial to generate the geometric phase for the multi-spin in a unit-cell unlike the one-spin case. It is 
found that the ground-state geometric phase, which is non-analytic at the critical points, possesses 
zigzagging behavior in the gapless B phase of non-Abelian anyon excitations, but is a smooth 
function in the gapped A phase. Furthermore, the finite-size scaling behavior of the non-analytic 
geometric phase along with its first- and second-order partial derivatives in the vicinity of critical 
points is shown to exhibit the universality. The divergent second-order derivative of geometric 
phase in the thermodynamic limit indicates the typical second-order phase transition and thus the 
topological quantum phase transition can be well described in terms of the geometric-phase. 

PACS numbers: 03.65.Vf, 75.10.Jm, 05.70.Jk 



I. INTRODUCTION 

Berry in his pioneer work raised a fundamentally im- 
portant concept known as geometric phase (GP) in ad- 
dition to the usual dynamic phase accumulated on the 
wave function of a quantum system, provided that the 
Hamiltonian varies with multi-parameters cyclically and 
adiabatically [l| . At the present time the GP with exten- 
sive generalization along many directions has wide appli- 
cations in various branches of physics ■ 

Recently, the close relation between GP and quantum 
phase transition (QPT) has been gradually revealed 0-01 
and increasing interest has been drawn to the role of GP 
in detecting QPT for various many-body systems [sl-fToj. 
which, as a matter of factj_js also a new research field in 
condensed matter physics • QPT usually describes 

an abrupt change in the ground state of a many-body sys- 
tem induced by quantum fluctuations. The phase tran- 
sition between ordered and disordered phases is accom- 
panied by symmetry breaking, which can also be charac- 
terized by Landau-type order parameters. 

On the other hand, a new type of QPT called topo- 
logical quantum phase transitions (TQPT) has attracted 
much attention. The first non-trivial example is the frac- 
tional quantum Hall effect [l^ [ij] . In the last decade, 
several exactly soluble spin-models with the TQPT, such 
as the toric-code model [lB|, the Wen-plaquette model 
16 , and the Kitaev model on a honeycomb lattice 
ISj , were found. In contrast to the conventional QPT 
governed by local order parameters the TQPT can 
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be characterized only by the topological order [19|. As 
good examples to illustrate the underlying physics, dif- 
ferent methods are developed to describe the TQ PT in 
the Kitaev honeycomb model [20l - [25j . In Ref. |2Cll |. Feng 
et.al. obtained the local order parameters of Landau 
type to characterize the phase transition by introducing 
Jordan- Wigner and spin-duality transformations into the 
Majorana representation of the honeycomb model. Gu 
et.al. showed an excitiiig result of the ground-state fi- 
delity susceptibility [l^, [1^ , which can be used to iden- 
tify the TQPT from the gapped A phase with Abelian 
anyon excitations to gapless B phase with non-Abelian 
anyon excitations. 

Quite recently, Zhu showed that the ground-state 
GP in the XY model is non-analytic with a diverged 
derivative with respect to the field strength at the critical 
value of magnetic field. Thereupon, the relation between 
the GP and the QPT is established. Nevertheless, much 
attention has been paid to the QPT, while effort devoted 
to the relation between the GP and the TQPT is very 
little. The present paper is devoted to exploiting the GP 
of the Kitaev honeycomb model as an essential tool to 
establish a relation between the GP and the TQPT and 
reveal the novel quantum criticality. Unlike the GP in 
the usual lattice-spin model for the QPT, which is gen- 
erated by the single-spin rotation of each lattice-site, the 
simultaneous rotation of linked two spins in one unit-cell 
seems crucial to describe the TQPT in the honeycomb 
model. The non-analyticity of GP at the critical points 
with a divergent second-order derivative with respect to 
the coupling parameters shows that the TQPT is the 
second-order transition and can be well described by the 
GP. 

In Sec. II, the ground state wave function and energy 
spectrum of the Kitaev honeycomb model are presented. 
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After introducing a correlated rotation of two z-link spins 
in each unit-cell, the ground-state GP and its derivatives 
are obtained explicitly in Sec. III. Sec. IV is devoted 
to investigating the scaling behavior of the GP. A brief 
summary and discussion are given in Sec. V. 



II. THE KITAEV HONEYCOMB MODEL AND 
SPECTRUM 

The Kitaev honeycomb model shown in Fig. [Ija) is 
firstly introduced to illustrate the topologically fault- 
tolerant quantum-information processing [ij, [l^ [2^ . In 
this model, each spin located at vertices of the lattice in- 
teracts with three nearest-neighbor spins through three 
types of bonds, depending on their directions. By using 
the Pauli operators cr" (a — x,y,z), the corresponding 
Hamiltonian is written as 



a;-links 



w- links 



J. 



2^ f^j^fc' 

-links 



(1) 



where j, fc denote the two ends of the corresponding 
bond, and Ja are coupling parameters. After introduc- 
ing a special notation Kjk ~ '^'j'^ki where the indexes a 
depend on the types of links between sites j and k (so 
we also write it into a^fe in the following text for per- 
spicuousness) , Hamiltonian ([1]) can be rewritten into a 
compact form 



U.k) 



(2) 




FIG. 1: (Color online) (a) Kitaev honeycomb model, in which 
one spin interacts with three nearest-neighbor spins througli 
three types of bonds, depending on their direction. A unit-cell 
with X, y and z links and graphic representation of Hamil- 
tonian Q with Majorana operators are marked by the red 
dotted line, (b) Phase diagram of the model, where A phase 
is gapped and B phase is gapless. 

It has been known that the Kitaev honeycomb model 
can be solved exactly by introducing Majorana fermion 
operators, which are defined as [3, 



ib^c. 



ib'c. 



(3) 



Generally, a set of Majorana operators M = {b^ ,by,b^,c} 
can be employed to describe a spin by two fermionic 



modes. They are Hermitian and obey the relations 
— 1 and mm' — —m'm for m,m' G M and m ^ m' . 
Moreover, in the Hilbert space with a spin described by 
two fermionic modes, the relation b^Wb^c \^) — |^) must 
be satisfied to ensure the obeying of the same algebraic 
relations as ct^, ct^, and ct^ |18| . 

Drawing on the operators to the Kitaev honeycomb 
model, Hamiltonian ^ is given by 



-y 



(4) 



where Uj^ — i6^''°fe^''°. Fig. [IJa) also shows the struc- 
ture of Hamiltonian (j4]), from which it can be seen that 
Ujk = —Ukj- Since these operators Ujk commute with the 
Hamiltonian (U) and with each other, the Hilbert space 
splits into two common eigenspaces of Ujk with eigenval- 
ues Ujk — ±1. Thus, Hamiltonian (j4]) is reduced to a 
quadratic Majorana fermionic Hamiltonian 



With a Fourier transformation 



1 iqr^ 



(5) 



(6) 



where s denotes a unit cell shown in Fig. [Ija), A refers 
to a position inside the cell, represents the coordinate 
of the unit cell, and q are momenta of the system with 
finite system-size 2L^, and a Bogoliubov transformation 



a 



'q,i 



q.2 - v5 



q,2, 



(7) 



where Aq = ^e^ + A'^/{Aq + ie^), Hamiltonian ([SJ is 
transformed into 

^ = E \/< + (c^Cqa - Cl.c^a) (8) 



with Cq = Jx COS Qx + Jy COS qy + Jz, and Aq — sin + 
Jy sin qy . In Hamiltonian ([8|) , the momenta take the val- 
ues [2411 



qx(y) 



2mT 



L - 1 



L - 1 



(9) 



when the system size is chosen as iV = 2L^ with L being 
an odd integer. Thus, the ground and the first-excited 
states are obtained by 



i*o>= n ^1.2 io) = n 7f (^'i^-q.i + io) , (10) 

q q ^ 

n <i = n 7f - ^>-q^2) io) , (11) 
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with the energy eigenvalues 

^o.i = ±E\/4 + ^q' (12) 
q 

It has been shown that the Kitaev honeycomb model 
([T]) has a rich phase diagram including a gapped phase 
with Abelian anyonic excitations (called A phase) and 
a gapless phase with non- Abelian anyonic excitations {B 
phase) [i^. In Fig. [TJb), the two phases A and B are sep- 
arated by three transition lines, i.e., Jx = 1/2, Jy — 1/2, 
and Jz = 1/2, which form a small triangle surrounding 
the B phase. Here, we only plot the energy spectrum 
([T2|) as a function of Jz for Jx = Jy (the vertical dot- 
and-dash line in Fig. [Ijb)) in Fig. [2] It can be seen 
from Fig. [5] that the energy-level degeneracy arises or 
lifts at certain points, which can be regarded as the pos- 
sible critical points of QPT [1, [HI. In Fig. [DJb, c, d), 
the degenerate points occur in the B phase, but disap- 
pear in the A phase as shown in Fig. [2jf). Moreover, the 
energy spectrum may have asymptotic degeneracy at the 
phase diagram edge seen from Fig. [2Ka) when the size 
of system tends to infinity. The non-analyticity points of 
ground state in B phase are actual level-crossing points. 

III. GEOMETRIC PHASE 

In general spin-chain systems characterized by 
Landau-type order parameters, the GPs can be generated 
by adiabatic and cyclic evolutions with a rotation of each 
spin coupled with its neighbors in the certain parameter 
spaces 6]. However, in the Kitaev honeycomb model, the 
two spins (Js^x and \ in the unit cell s is coupled with 
each other only by a z-link, as Fig. [Ija) shows. As a 
result, a simple rotation of a single-spin with the genera- 
tor (t| ^ can only " propagate" along the horizontal spin- 
chain through the x- and j/-link couplings but not along 
the vertical direction since the rotation-generator of a 
single-spin (t| _^ commutes with the z-link spin-operators 
in the Hamiltonian between the horizontal spin-chains. 
On other words, if using the above single-spin rotation, 
the Kitaev honeycomb model would be equivalent to a 
system of independent multi-single-spin-chains. In or- 
der to obtain the GP of the Kitaev honeycomb model, it 
should be introduced a correlated rotation of two z-link 
spins in each unit-cell, 

U (0) = exp [iipR] , (13) 

where 

^ = EE<A<A (14) 

is the correlated rotation-generator, and is the "co- 
rotation angle" . Indeed, when the rotation operator is 
applied on the spin of (s, A)-th set, the spin operator cr^ ^ 
acts as a c-number because the spin operators on different 




FIG. 2: Energy spectrum for the parameters ~ Jy (a) 
Jz = 0, the spectrum indeed is degenerate in a larger system 
size; (b), (c), and (d) in B phase, for Jz < 1/3, Jz = 1/3, 
and 1/3 < Jz < 1/2, respectively; (e) at the critical point 
of Jz — 1/2 and (f) in A phase. It is clearly that there is 
level-crossing in (b), (c) and (d). 

sets commute each other. As a matter of fact, the corre- 
lated rotation-generator R can generate a "co-rotation" 
of the z-link spins in each unit cell. Furthermore, under 
the gauge transformation ([Tn|) , the rotations of each spins 
in the 2D honeycomb Kitaev model become correlative 
and global. Hence it is possible to reveal the relation be- 
tween the GP and the nonlocal TQPT by means of the 
correlated rotation. Then, the ground state wave func- 
tion becomes 

!*[,(</.)) = C/(0)|*o), (15) 
and correspondingly, the GP is given by [27j 

Jo "0 
--27r(^'o|i?|^'o). (16) 

It can be seen that the GP for the Kitaev honeycomb 
model is proportional to the expectation value of the 
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correlated rotation R. Different from the local spin- 
correlation cr? i cr^ T , the sum over whole lattice-sites in 
R leads to a global property, which is crucial for the 
topological QPT. 

In terms of the Majorana fermion operators ([3]), the ro- 
tation operator R is given, after using the Fourier trans- 
formation dH), by i? = XIs.A {^K,\Cs,x) {^K.x^s,x) = 

Es (-«<l;s,2Cs,lC5,2) + (-«<2;s,lCs,2Ca,l). Under 

the restriction of vortex free subspace, Ug i.,, 2 u^^2-s.i 
are the good quantum numbers, that is to say, 2 — 
1, and 2 s 1 = Thus, the rotation operator R is 
finally obtained by 



R ^ {a^qAaq,-! - a^q,2aq,i) . (17) 

By means of Eqs. (fTO)) . (fT6)) . and p7)) . the ground-state 
GP is given formally by 



-27r 



(18) 



where d is a vector from one site to another inside one 
cell. It can be seen easily from Fig. [IJa) that the vector 
is given by d = (rii -I- n2) /3. Inserting this vector d into 
Eq. (dH), we have 



-27r 



^Im (-Aqet(29.+29«)+^; 



-27ry^ eq[l+cos|(2g ^+ 2^^)] -Aq sin^ {2q^+ 2qy) 



(19) 



L=101 




Eq. p9)) is the main result of this paper. The scaled 
(or average) GP 7/L as a function of coupling parame- 
ters is shown in Fig. |3l which distributes symmetrically 
with respect to the coupling parameters Jx and Jy. For 
sake of simplicity, the scaled GP 7/i is replaced by 7 in 
the following discussions. In order to see the difference 
between A and B phases, 7- variation with Jz along a se- 



and Jy, 
The 



FIG. 3: (Color online) Scaled GP with system size parameter 
L = 101. 



lected line Jx = Jy in the phase diagram of Kitaev [18 1 
(see Fig. [ijb), the vertical dot-and-dash line) is plotted 
in Fig. Ufa). It is shown that 7 is smooth in A phase 
(Jz > 1/2) while becomes saltant in B phase (J^ < 1/2) 
for the size parameters L = 11 (dark yellow line), 33 
(red), and 99 (blue), respectively. Moreover, all the data 
fall onto a single curve in A phase, while the number of 
saltation increases with the system-size L in B phase (see 
insets (1) and (2) of Fig. IUJa) ). To be specific, the num- 
ber of saltation for L = 33 is three times than that for 
L = 11 (see Fig. Ufa)), and the same situation occurs in 
turn for L = 99 and 33. 

It is meaningful to consider the first-order partial 
derivative = dj/dJp {/3 = x, y) of the GP 7. Since the 
GP 7 in Eq. (jl9|) is symmetric with respect to J; 
we only need to investigate gx (or equivalently gy) 
variation of gx with respect to Jz along the selected path 
of Jx = Jy from B phase to A phase is shown in Fig. 
|4jb) for different system-size parameters L — 101, 303 
and 909. It can be seen from Fig. IHb) that gx oscillates 
in B phase with frequency (or number of peaks), which is 
proportional to L (see inset (1) of Fig. (Hb)). This rapid 
variation of the GP 7 (in the B phase) has not yet been 
found, to our knowledge. However, a very similar behav- 
ior of fidelity susceptibility in the Kitaev model has been 
reported [2J]. On the other hand, the value of gx at the 
critical point Jz ^ 1/2 increases with the system-size and 
sharply decays in A phase (see inset (2) of Fig. IHb) for 
detail). It is interesting to remark that the saltation of 
GP 7 in i? phase due to the complex structure of ground 
state with degeneracy (Fig. [2] (b),(c),(d)) is not random 
rather has regulation, especially it tends to a regular os- 
cillation above the point, Jz = 1/3 , (see inset (1) of Fig. 
IHb)). The oscillation frequency depends linearly on the 
system size. 

To show the non-analyticity of GP at the critical 
points explicitly the second-order derivative gxx of the 
GP 7 with respect to the coupling parameters is calcu- 
lated. Fig. shows the variation of gxx with respect 
to Jz along the variation path of Jz — Jy for different 
system-size parameters L = 101, 303 and 909. Inset (1) 
reveals the increase of peak-number of gxx in B phase 
along with the system-size L similar to 7 and gx in be- 
havior. The second-order derivative gxx is divergent at 
the critical point Jz — 1/2 as shown in inset (2) indicat- 
ing that the TQPT of the Kitaev honeycomb model is a 
second-order transition, While the QPT of the XY spin 
chain Q and the Dicke model 0, Q has been shown to 
be the first-order transition with the divergent first-order 
derivative of the GP. We conclude that the non-analytic 
GP 7 can very well describe the TQPT in terms of the 
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J J J 

z z z 

FIG. 4: (Color online) (a) 7 curve along the selected variation path = Jy for L = 11, 33, and 99. Both insets reveal the 
increasing number of saltation in B phase proportional to the system size parameter L. (b) Qx and (c) gxx as a function of Jz 
along the variation path Jx = Jy for system size parameters L = 101, 303 and 909. The two insets are local enlarged-pictures, 
which show the vibration in B phase and the circumstances in the critical point respectively. 



Landau phase-transition theory. 

Similarly, we can also choose the variation path as Jz = 
1/4 (dashed line in Fig. [Ijb)) with two critical points 
Jx = 1/4 and Jx — 1/2. Qualitatively similar results are 
shown in Fig. [S] for different size parameters L = 101 
(red line), 303 (blue) and 707 (dark), respectively, where 
7-plot is a smooth curve in A phase for Jx < 1/4 or Jj; > 
1/2 and becomes saltant in B phase when 1/4 < J^; < 
1/2 (Fig. [Sfa)). The plots of gx and gxx are respectively 
shown in Fig. [Sfb) and (c) displaying the clear vibration 
in the B phase, g^ has a sharp peak at the critical points 
and gxx is divergent in the thermodynamic limit showing 
the characteristic of second-order phase transition. Inset 
of Fig. \5lc) shows the detail of the vibration in the B 
phase with the peak-number proportional to the system 
size. 
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FIG. 5: (Color online) (a) 7 curve, which is smooth in A phase 
and saltant in B phase, (b) gx curve and (c) gxx as a function 
of Jx for Jx varying along the path = 1/4, with critical 
points at Jx ~ 1/4 and Jx ~ 1/2. 



IV. FINITE-SIZE SCALING 

In order to quantify the non-analytic nature of GP 
through the critical points, and to further understand the 
quantum criticality, we investigate the scaling behavior 
of the GP and its first-order and second-order derivatives 
by the finite size scaling analysis j28]. Since the GP of 
a many-body system is an extensive quantity, it usually 
depends on the system size in the non-critical region. In 
the vicinity of critical point (J^ = 1/2 here), the GP 7 
follows a power law|29j 

7CX I J,- Jf I ""^ (20) 

with a being corresponding exponent. \Jz — J^\- 
dependence of 7 in the left-hand side of critical point ( 
Jz < Jz ) and right-hand side ( J2 > Jf^) is respectively 
plotted in Fig. El^a) and (b) for different system-size pa- 
rameters L = 101, 303 and 909. The corresponding ex- 
ponents a- = -0.99934 ± 0.00033 (left-hand side) and 
a+ = -0.83538 ± 0.00008 (right-hand side) are obtained 
from Fig. [SUa) and (b). Similarly, gx and gxx as a func- 
tion of jJ^ — J^\ are plotted in Fig. [D^c), (d) and (e), 
(f) for the gapless and gapped phases respectively, from 
which the critical exponents = —0.17529 ± 0.02056, 
a-^^ = 1.32523 ± 0.01719 (left)" for B phase and a+ = 
-0.48375 ±0.00033, = 0.60687 ± 0.00526 (right)"' for 
A phase are found. The fact of negative exponents a^, 
Qfg^ and positive a^^^ indicates that 7 and gx are finite 
while gxx is divergent at the critical point in the thermo- 
dynamic limit. Thus the TQPT is a second-order phase 
transition characterized by the GP 7. 

The position of maximum value of 7 denoted by Jip 
may not be located exactly at the critical point = 1/2, 
but tends to it in the thermodynamic limit L — cxd, which 
is regarded as the pseudocritical point [111 . The J^' — J™ 
versus different system-size L = 901, 951, 1901 in a 
logarithmic coordinate is plotted in Fig. [3 which is a 
straight line of slope —0.37851 ± 0.00226. It means that 



6 



(a) 




^-^^■""(2 L=901 
|> L=1201 
<] L=1501 
y=11.36197'|J^-J^|°^ 


^^-A. 26603 


(b) 










O L=901 ^^""^'^f-^s^ 
L=1201 
<] L=1601 

T=-3.70219*|J_-J°|°™''-4.26456 


- 











1 ° 

(c) 


o 

o 


O L=901 
L=1201 
<j L=1501 
g =-3.29946'| J^-J^l" 

^^"ft^'tM^ p ^ r r o o o c 


-21 .33245 








(d) 




Jt-^^'iS L=901 
|> L=1201 
<] L=1501 
g=100.54205*|J^-J°|™ 


^-22 3661 5 


9 

\ 




O L=901 
|> L=1201 
<\ L=1501 
9^=0.04982*1 J^-J°r' ™ 

-C . r - „ „ --J ^ 


-52.09722 


(e) 








(f) 




O L=901 
t> L=1201 
<\ L=1501 
g__=56.09834*|J_-J°r°" 












0.000 




0.002 0.004 


0-006 



FIG. 6: (Color online) Finite-size scaling analysis of the 
power-law divergence for (a) GP 7, (c) Qx, (e) g^x as a func- 
tion of I Jz — in the vicinity of critical point with system 
sizes L = 301,901 and 1501 on the left-hand side (J^ < jf ) 
and (b), (d), (f) on the right-hand side (J^ > J^) respectively. 

the tends toward to the critical point following 
the power-law decay 

Jf - jrocL'0-37«5i^ (21) 

On the other hand, the maximum value of 7 at = 
J™ for a finite-size system behaves as 

7(Jr)cxL^-, (22) 

which is shown in the inset of Fig. [SJa) with a straight 
line in logarithmic coordinate. The corresponding size- 
exponent is given by fij = 0.01148 ± 0.00054. 

Since the GP 7 around its maximum position J™ can 
be written as a simple function of J™ — Jz, it is possible 
to make all the value-data defined by a universal scal- 
ing function Fj = (7(J™) — 7) /"f versus L'^'' (Jz — J™), 




-5.30- 



6.8 6.9 7.0 7.1 7.2 7.3 7.4 7.5 7.6 

In L 

FIG. 7: (Color online) L-dependence of — J™ in logarith- 
mic coordinate for L = 901, 951, 1901. 

namely, 

F, = f[L-^^{JT-J.)], (23) 

where v-y is a critical exponent that governs the diver- 
gence of the correlation length. The values of for dif- 
ferent system-size parameters L fall onto a single curve 
as shown in Fig. ^sl), from which we can extract the 
critical exponent = —0.015 numerically. 

In fact, according to the scaling ansatz of a finite sys- 
tem [28l . [2^ , the critical exponent v can be determined 
by the relation — /i/a. In terms of this relation, 
the critical exponents in B and A phases are found as 
= —0.01149 and = —0.01374, which arc consis- 
tent with the numerical result Vj extracted from Fig. 
^a). Inset of Fig. [5]Jb) is a plot of the maximum value 
of g-xiJ"'') ^ function of L 

g^iJT) cx L'^- (24) 

in logarithmic coordinate, from which the size exponent 
fig^ = 0.01612 ± 0.00094 is found. The universal scahng 
function 

Fa.^f[L''-^{JT~Jz)] (25) 

is shown in Fig. [IJb). We have the numerical value 
j/g^ = —0.040 and the results determined by the relation 
1/ = fi/a that f- = -0.09196, j/+ = -0.03332. The de- 
viation between Vg and v^^^ may be due to the rapid os- 
cillation in the gapless B phase. Similarly, the critical ex- 
ponents of gxx can be obtained from Fig. [H^c) as fj,g^^ = 
0.54570 ± 0.00237, i/" = 0.410 and j/+ = 0.908, Ind 
the results determined hy v = ^/a are v^^^ — 0.41178 
and f ~ 0.89920 respectively. 

V. SUMMARY AND DISCUSSION 

We demonstrate that the ground-state GP generated 
by the correlated rotation of two linked-spins in a unit- 
cell indeed can be used to characterize the TQPT for the 
Kitaev honeycomb model. The non-analytic GP with a 
divergent second-order derivative at the critical points 



L^(J,-J") 



L^'<Jz-Jz) 



FIG. 8: (Color online) (a) F^, (b) Fg^, and (c) Fg^^ as a function of L" (J^ - JT) for L = 401, 601, 1601. All the data fall 
on to a single curve respectively. And the three insets in turn shows the variation of 7(J™), QxiJT), and gxx{JT) with respect 
to the system size-parameters L = 101, 151, 1901. 



shows that the TQPT is a second-order phase-transition 
different from the XY spin-chain ^6|, iSj], in which the 
first-order derivative of GP is divergent, and the LMG 
modeljl, Q, in which the GP itself is shown to be di- 
vergent. Moreover it is found that the GP is zigzagging 
with oscillating derivatives in the gapless B phase, but is 
a smooth function in the gapped A phase. The scaling 
behavior of the non-analytic GP in the vicinity of critical 
point is shown to exhibit the universality with negative 
exponents of both 7 and Qx while a positive exponent 
of Qxx indicating the characteristic of second-order phase 



transition. 
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